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Rossby wave turbulence, as modelled by the Charney-Hasegawa-Mima (CHM) equation, is nonlo- 
cal in scale. As a result, the formal stationary Kolmogorov-Zakharov solutions of the Rossby wave 
kinetic equation, which describe local cascades, are not valid. Rather the solution of the kinetic 
equation is dominated by interactions between the large and small scales. This suggests an alterna- 
tive analytic approach based on an expansion of the collision integral in a small parameter obtained 
from scale separation. This expansion approximates the integral collision operator in the kinetic 
equation by anisotropic diffusion operators in wavenumber space as first shown in a series of papers 
by Balk, Nazarenko and Zakharov in the early 1990's. In this note we summarise the foundations 
of this theory and provide the technical details which were absent from the original papers. 



I. INTRODUCTION TO WEAK WAVE 
TURBULENCE IN THE 
CHARNEY-HASEGAWA-MIMA EQUATION 

In the limit of weak nonlinearity, the statistical evolu- 
tion of CHM turbulence can be described by the theory 
of weak wave turbulence. See [1, 2] for a review of the 
theory. The principal result of this theory is the fact 
that that the wave spectrum, nk, evolves according to 
the wave kinetic equation: 



dt 



[nqTir — nknqSgn(cJkWr) — JT-k^ir sgn(a;ka;q)] dqdr. 
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Throughout these notes we shall use boldface text to de- 
note vectors in R^, regular text to denote their magni- 
tudes and subscripts to denote their components. Thus, 
for example, k = (fci,fc2) and /c = |k| = •\/fcf + AJf. The 
nonlinear interaction coefficient, has several symme- 
tries which will be useful later: 



V, 



(fel,fc2) 



(91,92) (ri,r2) 



^(ri.r^) (131,92) 

_T/(fel,-fe2) 

(91,-92) {ri,- 

y(fel,'£2) 
( — 91,92) (ri,r2) 



(ri,r2) 



(fel 



92) 



^2) (-91 

(91,92) (ri,r2) 
(kiM) 



(91 ,92) (-ri,r2) ' 



(4) 



•Electronic address 
t Electronic address 
t Electronic address 



coimaughtoncOgmail . com 

S . V . Nazarenko® . Warwick .ac.uk 

B .E. QuiimOwarwick. ac .uk 



FIG. 1: Typical shape of the resonant manifolds defined by 
Eq. (5). 



According to Eq. (1), exchange of energy is only possi- 
ble among triads of modes, (k, q, r), which satisfy the 
resonance conditions: 



k = q + r 

Wk = + UJr 



(5) 



For each given k, the modes q permitted to interact with 
k lie on a one-dimensional curve in the {ki,k2) plane 
known as the resonant manifold of the mode k (the third 
member of the triad is given by r = k — q) . The resonant 
manifolds typically have the shape shown in Fig. 1. 

The kinetic equation, Eq. (1), becomes scale invariant 
only in the limits, k ^ p^^ and k ^ p^^. In these lim- 
its, it is possible to formally find stationary Kolmogorov- 
Zakharov (KZ) solutions of Eq. (1) describing turbulent 
cascades of energy and enstrophy. The inverse cascade 
spectrum, describing the transfer of energy from small 
scales to large, turns out to be non-local meaning that 
the RHS of Eq. (1) diverges on this spectrum so that 
it is not a consistent stationary solution. It is worth 
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stressing that the fact that the "usual" KZ spectrum is 
not realisable docs not necessarily invalidate the use of 
the kinetic equation itself. Rather, the non-locality of 
the KZ solution leads one to expect that the transfer of 
energy from small scales to large in the CHM equation 
proceeds by a direct interaction between the large and 
small scales rather than via a scale-local cascade process. 
This is what is meant by non-local turbulence. In this 
section, we provide a detailed description of this process. 
This analysis was originally presented in [3, 4] and is re- 
produced here in slightly expanded form for the sake of 
clarity and completeness. 

The x-component of the phase speed of a wave having 
wavenumber k is 



Ck = 



Wk 

h 



(6) 



Notice that all waves have a phase velocity to the left. 
Furthermore for large scale waves, k Ig p~^, Ck — >■ —/3p^ = 
—V* and all waves travel at the same speed, known as the 
drift velocity: 



(7) 



We are interested in describing a scenario in which small 
scale turbulence evolves by direct interaction with large 
scale waves. Since the large scale waves are all travelling 
at the drift velocity in the x-direction, it is natural to 
work in a frame moving in the x-direction with the drift 
velocity, —(3p^. In this frame, using Eq. (6), the phase 
speed of mode k changes to: 



Ck 



Ck + /3p^ = 



The frequencies are Doppler-shifted accordingly: 
Wk = Ck fci ^ /3 ki + Wk- 



(8) 



We shall use flk to denote the frequency in the moving 
frame: 



Ok = P ki+LOu 



1 + p^k^ 



(9) 



In the moving frame the kinetic equation, Eq. (1), be- 
comes 



^ = 4TrJ\V^^f S{k-q-r)S{uJk-ujq-uJr)x 

[nqUr - TikTiq sgn(a;kWr) - nkTir sgn(wkWq)] dqdr, 

(10) 

Note from Eq. (8) that the change from Wk flk does 
not change the resonant manifolds, Eqs. (5), due to the 
presence of the momentum delta function, (5(k — ki — k2). 

Let us now consider small scale turbulence with typical 
wave- vector k. From Fig. 1, we can see that the resonant 
manifold containing the set of wave- vectors, ki, permit- 
ted to interact with k typically intersects the ki = axis 



in two places. Setting fci = in Eqs. (5) and doing some 
calculations shows that these points are 



Pi = 

P2 = 



(0,0) 

(0,-2A;2). 



(11) 
(12) 



If we are interested in non-local interaction between 
wavenumbers in the neighbourhood of k and zonal flows, 
then the points Pi and P2 will give the dominant contri- 
butions. Interaction with the point Pi always describes 
interaction with a large scale zonal flow. The wave num- 
ber of P2, on the other hand, is comparable to the fc2 com- 
ponent of the small scale turbulence. Interaction with 
the point P2 can correspond to interaction with a large 
or small scale zonal flow depending on whether the small 
scale turbulence has an appreciable ^2 component or not. 
We focus on the interaction with large scale zonal flows 
flrst. 



II. NONLOCAL INTERACTION WITH LARGE 
SCALE ZONAL FLOWS 

Let us introduce a large scale reference scale, K <^ p~^, 
and after integrating out r from Eq. (10), split the right 
hand side of the kinetic equation as 



dt 



+ 



I CoU[nk,nq,k, q] dq (13) 

Jq<K 

I Cofl [nk, riq, k, q] dq, 

Jn>K 



lq>K 

where the collision integral is: 

Coll [rik, nq, k, q] = 47r iK^k-qf ^(^^k - Oq - Ok-q) 

[riqrik-q - nknqSgn(a;kWk-q) - nknk-qSgn(wfeWq)] rfki 

(14) 

If we are interested in the dynamics of the small scales, 
k ^ K, then the assumption of nonlocality means that 
we can neglect the second term in Eq. (13) compared 
to the first. Further simplifications can be made in the 
first term using the fact that k ^ q everywhere in the 
integrand and assuming that Uk <C riq everywhere in 
the integrand. This latter inequality means that we are 
describing a situation where small scale turbulence is in- 
teracting with an intense large scale zonal flow. This 
discussion does not tell us how this large scale zonal flow 
was generated in the flrst place. 

As q -> 0, sgn(cLik'^k-q) = 1 and wc can approximate 
the rik dependence in the coUison integral as follow: 

[nqTT-k-q - ^^k^^qSgn(wkWk-q) - nknk-qSgn(wfeWq)] 

-> ("k-q - «k) nq 

where we neglect the term nk?^k-q describing small-scale 
small-scale interactions on the basis that it is O(n^) as 
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q — )■ and much smaller than the terms describing small- 
scale large-scale interactions. In other words, we assume 
that rik <^ riq. We therefore approximate Eq. (10) by 



dt 



= I i^(k,q) 

Jq<K 



rfq 



(15) 



where 



J^(k,q) 



T/1 



47r 

xnq (nk-q 



(5(r2k - - ilk-q) 



(16) 



Using the symmetries of Vj^^ and Ok it can be shown 
by direct computation that 



i^(k,q) = 
Thus we can then write: 



-i^(k-q,q). 



(17) 



dt 



rfq 



= \ ( F(k,q)rfq- / F(k-q,-q; 

^ Jq<K Jq<K 

= \ j (F(k,q)-F(k + q,q)) dq 

^ Jq<K 

= -\ I dqq-VkF(k,q), (18) 

^ Jq<K 

where in the final step wc have Taylor expanded _F(k + 
q, q) with respect to q in the first argument and neglected 
terms of 0{q^). 

Wc shoiild now Taylor expand -F(k, q) with respect 
to q in Eq. (18) To do so, we first Taylor expand the 
argument of the ^-function. Taylor expanding Ok-q with 
respect to q gives: 



O] 



k-q 



Ok - q • VkOk + 0{q^ 



as q 0, where Vk = {dk^ , ) is the gradient operator 
in k-space. Thus the difference Ok — Ok-q behaves as 
q • VkOk as q — > 0. On the other hand, the remaining 
term Oq behaves as q^qi as q — >■ as can be seen directly 
from Eq. (9). Thus the latter can be neglected compared 
with the former in the limit q ^ and we can replace 

5(0k - Oq - Ok-q) -)> 5(q • VkOk) 

in Eq. (14). Therefore F(k, q) in Eq. (18) can be approx- 
imated as 

F(k,q) w 47r|Vq''k-qr'^(q- VkOk)(nk-q-nk)nq 

~ -47r|yq\_qf (5(q-VkOk)q-Vknk, (19) 

where we have used the fact that rik-q — ?^k = — q • 
Vk^k -I- O(g^) as q — 0. Combining Eq. (19) with 
Eq. (18), the kinetic equation for the small scales can be 
written as an anisotropic diffusion equation in k-space: 



d 



(20) 



where repeated component subscripts are summed over. 
The diffusion tensor is 



5i,(k) 



27r 



dq 



q<K 



q k-q I 



(5(q-VkOk)g'igjnq. (21) 



Let us now suppose that the large scales are principally 
supported at scales q <^ K. The reference wavenumber, 
K, can then be extended to infinity in Eq. (21) and we 
can perform the integration with respect to qi : 



o I J J iT^k |2./ 90k , 90k, 

27r / ctgia«2 |Kj k-q| '^('^'i 5^ +«2-^)Q'i9jnq 



Sij{]f.) = 2n f 

J — C 

/oo 
-OO 



2n 



90k 



dki 



f 

J — ( 



dq2 



90l 



dki 



S{qi +0kq2)qiqjnq 



Ik" 



qk-q| QiQjn^ 



qi= — 6k«2 



where we have introduced the shorthand parameter 

- on^ 

dki 

Some elementary calculations bring the diffusion tensor 
to the form: 



(22) 



5(k) = 5(k) 



90k 



01 
-9w 



(23) 



where we have introduced the scalar function 



5(k) = 2tt 



dq2 



I "^q k-q| 



91 = 



(24) 

The matrix in Eq. (23) is singular. It has eigenvalues 
and 1 + ^k ^ith corresponding eigenvectors (^k^, 1) and 

{—0, 1). This suggests that at any point in k-space, diffu- 
sion only occurs in the direction (— 0k, 1) corresponding 
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to the positive eigenvalue. Eq. (20) should then describe 
diffusion along one-dimensional curves in k-spacc. In or- 
der to show this and to find these curves, we need to 
find a change of variables which reduces Eq. (20) to a 
1-dimensional diffusion equation. 
The general structure of Eq. (20) is 

^=(^.^)^«(t)- <-) 



Under a change of coordinates, 

fci ^ Kl(fci,fc2) (26) 
k2 K2{ki,k2), 

Eq. (25) transforms to [5] 



dt 



Idet J| 



J5(k) J' 




where our notation tacitly assumes that Eq. (26) is used 
to express all k dependence in S'(k) and rik in terms of 
Ki and K2- J is the Jacobian matrix of the change of 
variables: 



(28) 



For a typical anisotropic diffusion equation for which the 

diffusion tensor would have two positive eigenvahies, the 
maximum simplification would be obtained if new coor- 
dinates could be found so that J ^(k) is diagonal. In 
this case, since the diffusion tensor is singular, we can 
do better and reduce the equation to a one-dimensional 
diffusion equation. By direct computation we observe 
that 



a b 
c d 



01 



-Ov 1 ) \ c d 

{h-aOi^f {b-a 
{b-a9i,){d-c9i,) 



Using this we see that 



01 



-Ok 

1 



a9]f^){d — ct 
{d-c9i,f 



1 




provided that our new coordinates satisfy the equations 

dni _ dni _ ^ 
dk2 dki 

dK2 _ dK2 ^ Q 

dk2 dki 

The first of these equations can be easily solved by in- 
spection. The second is also obvious once one recalls the 
definition of 9\^ from Eq. (22). We obtain: 



d 
a 

dK,2 



(27) 



r 




FIG. 2: 

The level sets of the function fik given by Eq. (9). 

For this change of variables, det J = so that in the 
(ki,k2) plane, the diffusion equation Eq. (25) takes the 
form 



dki dKi Oki 



(30) 



In the (k;i,K2) plane Eq. (30) describes one dimensional 
diffusion in the ni direction with K2 constant. Using 
Eqs. (29) to translate this back into the (fci,fc2) plane, 
Eq. (25) therefore describes one- dimensional diffusion in 
the k2 direction along lines of constant fik- Illustra- 
tive examples of the curves Ok =constant are plotted 
in Fig. 2. 



III. NONLOCAL INTERACTION WITH SMALL 
SCALE ZONAL FLOWS 



Kl(fci,fc2) = k2 

K2{ki,k2) = 0(^1,^2). 



(29) 



In this section we perform the analogous calculation 
assuming that the evolution is dominated by interaction 
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with the point P2 = (0, 2fc2) corresponding to smaU scale 
zonal flows. The results of this section were first pre- 
sented in [6]. 

Following the same reasoning as before, we introduce 
a reference scale, K, satisfying K <^ k and isT <C 1/p, 
and neglect the contributions to the collision integral 
coming from outside of a ball of radius K of the point 
p = (0, 2 fc2)- We thus approximate Eq. (10) by 



dt ~ Ju 



|p-q|<if 



dqF(k,q) 



|p-q|<K 



dqF(k-q,-q). (31) 



Taylor expanding the first argument about q = p gives 



9nk _ _ f 



dq 



|p-q|<K 



F(k - p, -q) + (q - p) • Vq, i^(k - qi, q)!^^^^ + 0(|p - q]') 



dq 



'|p-q|<K 

where we have introduced the shorthand notation 



F(k - p, -q) - (gi, 92 - 2 fe) • (afe,F(k, -q), afe,F(k, -q)) 



(32) 



k = k-p = {ki,-k2). 



(33) 



We now expand ^(k, — q) about q = p. Since Iqk-qk and riq both vary rapidly near q = p, we should only expand 
the ?T-i5._|_q — JT-k and 6{fl^^ — fl-q — ^^k+q) t^rms. Let us look at these two terms in turn. 



'^k+q-'^k = '^k+p-'^k + (q-P)-'^q' "(k + q') , +C'(|p-q|^) 

q'=p 

~ rik - "k + (91, 92-2 k2) ■ {dkirik, Sfe^rik) • 



(34) 



Note that, unlike the expansion of the collision integral about the point (0, 0) corresponding to large scale zonal flows 
detailed in Sec. II, the leading order term in the Taylor expansion of F{h, — q) about q = p is not necessarily zero. 
It vanishes only if the spectrum is symmetric about the ki axis (ie n{k\,k2) = n{ki, — ^2)). We shall return to this 
point below. Next, let us look at the argument of the 5- function near q = p. 



— fl-q — f^k+q ~ ^k ~ ^-P ~ + Ql 



-92) + ^(^1 + -^2 + g^) 



(9;,<if^)=(0,2fc2) 



+(g2-2A;2) 



+ 0(|p-q|') 



(«l,g2)=(o.2fc2) 



dk2 



(g2-2A;2+^kgi), 



where 



3(fcf - kj) + p^ikt + Qklki - 3A;|) 



2fcifc2(l+Vfc2) 

The various derivatives have been computed from Eq. (9). For example: 



(35) 



(36) 



7iki+q[,-k2 + q'2) 



W,92) = (0,2fe2) 



'dk^ 



2/3/fcifc2 

(l+p2fc2)2' 



and so forth. The leading order term is zero owing to the fact that k, p and k are resonant. Putting this together, 
the Taylor expansion of the (5-function is 



5(nfc - n_q - Ofc+q) 



dk2 



5 {q2 - 2A;2 + ^k9i) • 



(37) 
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From Eqs. (32), (34) and (37) we see that the leading order term in the kinetic equation coming from nonlocal 
interaction with small scale zonal flows is: 



= Yy, [n{ki, -k2) - n(fci, ^2)] , 



(38) 



where 



477 



= 477 



dk2 



dqidq2 



p-q.\<K 



V 



qk+qk 



dk2 



dko 



I 

J\p- 
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-<l\<K 



6{q2 -2k2 + ^kqi)n{qi,q2) 
dqidq2 |- Kik-^ikT - 2 ^2 - ^ qi) ■n{-qi, ^2) 
dqi. 



|V"qk-qk| nq 



92=2 k-2 



(39) 



In the intermediate step, we have used the symmetries, Eq. (4), of V^^^ and relabelled the integration variable 
q\ — >■ — gi to bring the interaction coefficient to a neater form. At the final step the delta function has been used to 
integrate out 52 to leading order. We have also assumed that the spectrum is symmetric about the k2 axis although 
this is an inessential point. Eq. (38) tells us that the leading order effect of interactions with small scale zonal flows is 
to cause the spectrum to relax to a state which is symmetric about the k\_ axis, a point which was first made in [3, 4]. 

If we wish to observe any redistribution of spectral energy density due to the interactions with small scale zonal 
flows, we need to consider the higher order terms in Eq. (32). The first order term vanishes after integration over q 
since the integrand is an odd function of p— q. The next contribution is therefore the second order one. From Eqs. (32) 
and (34) the second order contribution can again be presented as an anisotropic diffusion equation in k-space: 



9nk 



dt dki 



^ D n ^ ^'^k 
Bij(k)- 



dkn 



(40) 



where the diffusion tensor is given by the matrix 



B{\l) = An 



47r 



dk2 



dqidq2 



dko 



\p-q\<K 

/ dqi 

J|p-q|<K 



V 



qk+qk 



5{q2 - 2 fc2 + ^k qi) n{qi, 92) 



qi{q2 - 2 fe) 
91(92-2^2) {q2-2k2f 



dq2 l-V^jk-qkl 5{q2 - 2 ^2 - ^ ^i) «(-ai, ^2) 



-qi{q2 - 2^2) 
-qi{q2-2k2) (?2-2fc2)2 



= 47r 
9flk 



afik 


-\ 


1 -CkA 




dk2 




-a ) 


J [ 

J\qi\<K L 



|Vqk-Kjk| n{-qi,q2)ql 



92=2 k2 



dqi 



dk2 



B{k) 



-a 



(41) 



where we have performed the same manipulations as those used to arrive at Eq. (39) above and defined the scalar 
quantity 



B(k) 



= 47r / 



\<K 



|Vqk-qk| n{-qi,q2)ql 



g2=2/s2 



dqi. 



(42) 



The diffusion tensor, B(k), is again singular indicating that Eq. (40) should be reducible to a one-dimensional diffusion 
equation by an appropriate change of variables. Following the same proceedure as in Sec. II, we introduce new variables 



ki 
k2 



Ki{ki,k2) 

K2{kl,k2), 



(43) 



whose Jacobian, J, given by 



I C)K2 OK2 
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should diagonalise the matrix 



a b 
c d 



we see that 



appearing in the diffusion tensor, Eq. (41). Observing that 



a b 
c d 



(a-b^k)^ 
(a-6^k)(c-d^k) 



(a-Kk)(c-rfek) 
(c-da)' 



J 



1 



provided that our new coordinates satisfy the equations 



r = 



1 




dki dk2 
dki dk2 



1 
0. 



The first equation can be easily solved by inspection while the second requires a little more effort using the method 
of characteristics. One obtains 



Ki(k) = ki 
K2(k) = = arctan 



^2 + ki 
pk"^ 



arctan 



^2 - \/3 fci 
p k^ 



_ 2y/3pki 

~ 1 + p2 fc2 • 

With this done, det J = so that, according to Eq. (27), the diffusion equation, Eq. (40), transforms into 



driK, dZ]^ d 
dt dk2 dK\ 



dz^ 

dko 



'dko 



(44) 



(45) 



(46) 



where B(k) is given by Eq. (42). In the (ki, K2) plane Eq. (46) describes one dimensional diffusion in the ki direction 
with K2 constant. Translating back into the (^1,^2) plane, using Eqs. (44) and (45), Eq. (46) therefore describes one- 
dimensional diffusion in the ki direction along lines of constant Zk. Illustrative examples of the curves Zk =constant 
are plotted in Fig. 3. 



Interestingly Zk defined in Eq. (45), is closely related 

to the third invariant (in additional to the energy and 
enstrophy) of the full kinetic equation, Eq. (1), discovered 
in [7] (see also [8, 9]). From Fig. 3 we note that the 
diffusion due to interaction with small scale zonal flows 
tends to transfer the wave action to large scales, a process 
which was investigated numerically in [10]. 

IV. CONCLUSIONS 

In this note we have summarised the derivation of 
the nonlocal kinetic equation for Rossby wave turbulence 
modelled by the wave kinetic equation obtained from the 
Charney-Hasegawa-Mima equation. For small scale tur- 
bulence there are two possible sources of nonlocal inter- 
actions in k-space. The first is with large scale zonal 



flows with spectral support around the point Pi = (0, 0). 

The second is with small scale zonal flows with spectral 
support around the point P2 = (0,2^2). In both cases, 
the redistribution of small-scale wave action in k-space 
is described at leading order by an anisotropic diffusion 
equation in k-space. In both cases, the diffusion tensor 
is singular meaning that, via an appropriate change of 
variables, it can be shown that the diffusion occurs along 
one- dimensional curves in k-space. We showed how to 
obtain these curves explicitly. In the case of nonlocal 
interaction with the point Pi, the curves are open and 
tend to transport waveaction to ^2 = ±00. In the case 
of nonlocal interaction with the point P2, the curves are 
closed and tend to transport wave action towards the ori- 
gin. Recent numerical investigations have demonstrated 
the presence of these diffusive mechanisms in simulations 
of the original CHM equation [11]. 
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FIG. 3: 

The level sets of the function given by Eq. (45). 
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